PERIODIC DE BRUIJN TRIANGLES: EXACT AND ASYMPTOTIC RESULTS 



B. Shapiro*, M. Shapiro* and A. Vainshtein^' 

^ Department of Mathematics, University of Stockholm 
S-10691, Sweden, shapiro@inateinatik.su.se 
* Department of Mathematics, Michigan State University 
East Lansing, MI 48824, USA, mshapiro@math.msu.edu 
^ Dept. of Mathematics and Computer Science, University of Haifa 
Mount Carmel, 31905 Haifa, Israel, alek@mathcsll.haifa.ac.il 

Abstract. We study the distribution of the number of permutations with a given periodic up-down 
sequence w.r.t. the last entry, find exponential generating functions and prove asymptotic formulas for 
this distribution. 



§1. Introduction and results 

Let a — ((Ti, . . . ,cr„) be a permutation of length n. We associate with a its up-down sequence 
(sometimes called the shape of a, or the signature of a) V{<7) — {pi, . . . ,p„-i), which is a binary 
vector of length n~ I such that = 1 if ct^ < Ci+i and Pi ^ otherwise. During the last 120 years, 
many authors have studied the number (t^ of all permutations of length n with a given up-down 
sequence V. Apparently, for the first time this problem was investigated by D. Andre [Anl, An2], 
who considered the so-called alternating (or up-down) permutations corresponding to the sequence 
7-* = (1, 0, 1, 0, . . . ) = (10)* and proved that the exponential generating function for the number of 
such permutations is equal to tan a; sec a;. In [An3] he proved that this number grows asymptotically 
as 2n!(2/7r)"+i. 

A general approach to this problem was suggested by MacMahon (see [MM]). This approach leads 
to determinantal formulas for jj^, rediscovered later by Niven [Ni] from very basic combinatorial 
considerations. For the relations of this approach to the representation theory of the symmetric 
group, and for its generalizations, see [Fo, Stl, BW]. 

Another, purely combinatorial approach to the same problem was suggested by Carlitz [Cal]. His 
general recursive formula for (t^ is rather difficult to use. However, he managed to obtain explicit 
expressions for the corresponding exponential generating functions for certain periodic cases, that 
is for up-down sequences of the form V = (p)*, where p is a binary vector of a fixed length called 
the period of V. In [Cal, Ca2] he considered the case V = (1*^0)* and expressed the corresponding 
generating function via the Olivier functions of the kth order 



Another case, V = (1^0^)* , was considered in [CS1,CS2] and solved via Olivier functions of the fourth 
order. It follows that asymptotically jj^ in this case grows as 4n!(2/7)"+^, where 7 = 3.7502 ... is 
the smallest positive solution of the equation cost cosh t -1-1 = 0. 

The general periodic problem was solved completely in [CGJN]. As in the two particular cases 
mentioned above, the answer is expressed via Olivier functions. The techniques used involves matrix 
Riccati equations, and is rather complicated. For a different solution based on Mobius functions see 
[St2, Ch. 3.16, and Ex. 3.80]. 
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An additional dimension in the problem was introduced by Entringcr [En] who studied the dis- 
tribution of the alternating permutations by the last entry. He observed that the number Uj^j of 
alternating permutations of length i whose last entry equals j satisfy the following equations: 

llij = tli,j-i + tti-ij-i; tli,i = 0, i = 2k,k>0, 
= +tJ^-i J, k^=0^ i = 2k + l,k>0, 

with Hi^i = 1. These equations can be represented graphically as the following triangle 

1 

1 
1 1 
12 2 
5 5 4 2 

Fig. 1. The Entringer triangle 

Each even row of the triangle starts with 0, and an entry in such a row is equal to the sum of 
its left neighbors in the current and in the previous rows. Similarly, each odd row (except for the 
first one) ends with 0, and an entry in such a row is equal to the sum of its right neighbors in the 
current and in the previous rows. 

The Entringer triangle was studied by many authors. In particular, Arnold [Arl, Ar2] gave an 
interpretation of the entries of this triangle in terms of real polynomials with real critical values. 
Besides, he considered the exponential generating function 

and proved that A{x, y) = / cosh(a; + y). In fact, A{x, y) is the generating function of the signed 
Entringer triangle, which is obtained from the ordinary one by reversing signs in each ith row, 
where i equals or 3 modulo 4. Observe that the entries %j of the signed Entringer triangle satisfy 
relations 

(1-2) = + 

with boundary conditions fli^i = for i = 2fc, fli^i = for i = 2fc+l, fc > 0, = 1. General triangles 
satisfying relation (1.2) with arbitrary boundary conditions were first studied more than 120 years 
ago by Seidel [Se]. In particular, he proved that the ratio of exponential generating functions for the 
numbers on the right and on the left sides of such a triangle equals e^. More recently such triangles 
where studied, from the combinatorial point of view, in [DV, Dul, Du2]. In particular, it is proved 
in [DV] that the exponential generating function for a Seidel triangle is equal to e'^F{x + y), where 
F{x) is the corresponding function for the left side of the triangle. 

The case of general up-down sequences was addressed by de Bruijn [dB] (see also [Vi] for another 
version of the same result). Let ^fj be the number of permutations of length i whose last entry 
equals j and whose up-down sequence equals V = {pi,p2, ■ ■ ■)■ He proved that these numbers satisfy 
the following equations: 

C,- =Ci-i+Ci,,-i. Ci=0' ifft-i = l 



with fill = 1. Evidently, for P = (10)* one gets the Entringer relations (1.1). As before, these 
equations can be represented graphically as a triangle, and the direction in which one has to advance 
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along the rows of the triangle is governed by the sequence V. We call this triangle the de Bruijn 
triangle corresponding to the up-down sequence 7-*. A de Bruijn triangle is said to be periodic if the 
corresponding up-down sequence is periodic. 

Let V be a periodic up-down sequence with period p of length m > 1, and let ii < ^2 < • • • < V 
be the locations of zeros in p. Without loss of generality we assume that ir = rn (otherwise we 
consider instead of V the up-down sequence V = (pi)* , where pj = 1 — for 1 < i ^ m; evidently, 
the de Bruijn triangle for V is obtained from that for V by the reflection in the vertical axis). 

The signed de Bruijn triangle is obtained from the ordinary de Bruijn triangle by multiplying its 
ith row by 



(1.3) 



i > 1. 



The corresponding exponential generating function is defined by 



(1.4) 



i^l 3 = 1 



Theorem 1. The exponential generating function of the signed periodic de Bruijn triangle corre- 
sponding to the up-down sequence V is given by F^{x, y) = e^ f'^{x -\- y), where 



det M'P(t) 
detM^(t) 



and M'^{t) and M^{t) are r x r matrices 



M^{t) = 



( fmfi 



Vmfi 
Vm,i2— n 



^m,m+ii—i2 
^m,ir-i—i2 



^m,m-t-ii— v-i 
^m,m-t-i2— ir-i 

<Pm,0 



and 



M^{t) = 



fm,j{t)- 



( 1 



1 1 

y^m,0 ^m,m-\-ii—i2 
^m,i2—ii V^m,0 

Pm,ir-i—ii ^m,i^ — i —12 



^m,m-\-i2 —ir-i 
^m,0 



with ifm 

In particular, for signed Entringer numbers one has m = 2^ r — 1, and hence /{i">* {t) = Lp^l{t) = 
1/coshf, thus recovering the Arnold formula for A{x,y). Moreover, the same techniques allows to 
obtain generating functions for other Seidel triangles with periodic boundary conditions, such as the 
triangle for Genocchi numbers (see [DV]). It can be also extended to pairs of Seidel triangles with 
periodic boundary conditions, such as Arnold triangles L{(3) and R{f3) for Springer numbers (see 
[Ar2, Du2]), thus recovering several combinatorial results obtained in [Sp, Ar2]; see §2 for details. 

As an immediate corollary of Theorem 1 we get the above mentioned results of [CGJN] and [St2] 
concerning the generating functions for the number of permutations with a given up-down sequence. 

Corollary 1. Let r be even, then the exponential generating function for the numbers jj^ is equal 



1 + 



det M'Pjt) 
detM^it)' 
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where M'^(t) is an r x r matrix 



i>m,l i>m,2 
fm,i2-ii Vrn,0 



with 



1pm,r-l \ 

^m.ra+ii — 7',._ i 

711 + 12— ir-l 

^m,0 



and (fimj = Vm,j(t). 

To get a similar result for odd r, it suffices to consider P as a periodic sequence with a period of 
length 2m. 

Let us now consider the asymptotic behavior of the numbers ^fj. It was observed without a 
proof in [Ar2, p. 18] that the even rows of the ordinary Entringer triangle approximate, after an 
appropriate normalization, the function sin a; on the interval [0, 7r/2], while the odd rows approximate 
cos a;. Exact statements with the first two correction terms can be found in [SY]. 

We generalize this result to arbitrary periodic de Bruijn triangles. 

Theorem 2. For any I, ^ I ^ m — 1, one has 



lim 



aV \mn+l 

{mn + l- 1)! 



where Cm,i is a constant depending only on m and I, and uf is the normalized first eigenfunction of 
the two-point spectral problem 

with m homogeneous boundary conditions 



«W(0) = 
= 



if Pm+l- 
if Pm+l- 



-1 = 1, 
-1=0, 



where r is the number of zeros in the period p. 

The first eigenvalue of the above spectral problem is the minimal absolute value among the solu- 
tions of the equation 

detM^(i) = 0, 

and the eigenfunction uf is normalized by the condition vnaxQ^t^iuf {t) = 1. 
In particular, for the Entringer numbers tii,j one gets the sine law of [SY]: 

, 2fc-|-l 



lim 
lim 



hk+l,j (f ) 

(2fc)! 



■Kt 

= C2,l COS 



.t {2k-l)\ 



= C2,o sm 



TTt 

2"' 



it is shown in [SY] that C2,o = C2,i = 2. 

The starting point of this research was the result by the first and the third author that the 
numbers for V = (1^0^)* arise naturally in counting real rational functions of a certain type, see 
[SV]. The authors are grateful to Max-Planck-Institut fiir Mathematik, Bonn for its hospitality in 
September 2000 and to the Royal Institute of Technology, Stockholm for the financial support of the 
visit of A. V. to Stockholm in July- August 2001. Sincere thanks goes to R. Ehrenborg, A. Laptev, 
H. Shapiro, P. Yuditski and A.Volberg for useful discussions. 
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§2. Proofs 

Proof of Theorem 1. Define j^^- = £ifi^j- Since by (1.3), equals ei-\ if = 1 and — e^-i 
otherwise, wc immediately get that jj^^ = ^fj^i + j-i- In terms of the generating fimction 
defined by (1.4) this relation translates to = dF^ jdy — dF^ /dx. The general solution of this 
equation is given by F^{x,y) = e^ f^{x + y), where is a function of one variable to be defined 
from the boundary conditions. Let FJ^ and F^ be the restrictions of i^^ to the left and the right 
side of the signed de Bruijn triangle. It follows from the above discussion that F^{t) = f^{t) and 
^r{^) — 6*/^(*)- Evidently, there exist unique representations 



i=0 



besides, e* — X]i=o^ 'Pm,i{t)- Combining these expressions together we get 



(2.1) 



/ 9^ in \ 

tgTin 



m—1 



I ¥'m,0 Vm,m-1 '^m,m-2 
V'm,2 Vm,! fmfi 



^m,m— 1 ^m.rn — 2 ^m.m — 3 



Vm,3 



tfTitn 
tiffin 



where ipm,i = Vm,i{t)- Besides, relations pi = for i = ii, . . . , v- 
fc = 0, 1, . . . , and hence gf = for i = ii, . . . , V-i. Similarly, pm 
9o = 1> and = 1 for i ^ ii, . . . , v imply = for fc = 0, 1, . 

i 7^ 0, ii, . . . , ir-i- Therefore (2.1) is reduced to 



^m,o' \t"'-'fZ-iin' 

1 imply tiL+i+i,fem+i+i = for 
= together with ^f-^ = 1 imply 





^ ^m,0 






^m,m—ir-i 


\ 









y^m,m-\-ii—i2 


V^m,m+ii — Zt-- 


1 


= 


- ^m,i2 






y^m,m-\-i2—ir- 


1 








^m,ir-i—i2 


^m,0 





and hence f^ = for 



t'^fZin 



and the result follows. □ 



The Seidel triangle {"fij} for signed Genocchi numbers presented in [DV] is satisfies periodic 
boundary conditions 72fe,2fe = for fc = 1, 2, . . . , 72,2 = 1 and 72^+1,1 + 72fe+i,2fe+i = for /c = 
0, 1, Thus we arrive to equations 



t 



V>2,l{t) <fi2,o{t) 



tflif) 



which mean that the generating function for this triangle is 2{x + ?/)e^/(e^"'"^ + 1). 

In a slightly different way one can treat signed versions of Arnold's pairs of triangles {i(/9), i?(/9)}, 
{L{b),R(h)}, and {L{d),R{d)} involving Euler and Springer numbers (see [Ar2, Du2]). The first 
pair consists of triangles and satisfying periodic boundary conditions 



'2fc+l,l 



0, 



1,2,..., 



fc — 0, 1, . 
I^2k,2k = 0> A; = 1, 2, . 
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Their signed versions obtained by multiplying the ith row by (— l)(*~i)(*~2)/2 ^j.g ggi(jgi triangles. 

Similarly, the second pair consists of triangles {bfj} and {b^j} satisfying periodic boundary 
conditions 

&2fe+l,l = 0, fc = 0,l,..., 

^k,k — "k,l' ft — Z, O, . . . , i — 1, 

^2fe,2fc =0) k = 1,2, . . . . 

To get a pair of Scidcl triangles one has to multiply the ?'th row by (—1)*'^*"-'^)/^. 

Finally, the third pair consists of triangles {dfj} and {d^j} satisfying periodic boundary condi- 
tions 

rf2fe,i=0, A- = 1,2,..., 

dk,k = c^M' k = 2,3, . . . , 

'^2fe+l,2fe+l — 'J) k = 1,2, . . . , df^i — 1. 

To get a pair of Seidel triangles one has to multiply the ith row by (— l)(*^i)(*^2)/2 

Let f(t) = f^if) + tfl(f) and g<\t) = gl{t^) + tg{{i^) be the restrictions of the generating 
functions for the signed versions of and -R(/3) to the left and right sides, respectively. Then we 
get the following equations: 

J-Us.iW V2,oWy' V*/iV)v'' 
and hence the generating functions for the triangles and -R(/3) are 

^^^""'^^^ cosh2(x + ?/)' ^fl(^'y) = cosh2(ar + y)- 
For similar restrictions /''(t) and g^(t) one has 



; V^2,iW ¥'2,ow; v*A'(*')y v^"'sinht^ 

and hence the generating functions for the for the signed versions of triangles and are 

e^sinh(a; + j/) e"" cosh(a; + y) 

cosh2(a; + y)' cosh2(x + y) ' 

Finally, for f^if) and ^'^(t) one has 

1 \ ^ /^2,o(t) V's.iWV , (^cosht 

tgi^e)) \ y^^sinht 

and hence the generating functions for the for the signed versions of triangles L{d) and are 

, _ eni-cosh(x + y)) , _ e^+^^(l-cosh(x + y)) 

" cosh2(a; + y) ' ~ cosh2(a; + 2/) " 

Proof of Corollary 1. It is enough to observe that for r even, the exponential generating function 
for tt^ is equal to 1 + X^^q^ Si+it^ifT + sD' ^''^'^ to evaluate this sum according to the proof of 
Theorem 1. □ 
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Proof of Theorem 2. Consider the space PC[0, 1] of piecewise continuous functions on the segment 
[0,1] with the norm ||/|| = maxo^t^i Define two operators, S° and S^, taking PC[0, 1] to 

itself: 

S'fix) = [ f{t)dt, S'fix) = [ f{t)dt. 







It is easy to see that both S° and can be written as integral operators 

S'f{x)= f K\x,y)f{y)dy, z = 0,l, 
Jo 



with the kernels 

K'^ix.y) 



1, x^y, 
0, X <y, 



and K^{x, y) — I — K'^{x, y). Besides, consider two families of integral operators S*^, i = 0, 1, with 
the kernels K'^{x, y) given by 

A„(a;, y) = s 

y 0, otherwise, 

and K^{x,y) = 1 — K^{x,y). It is easy to see that 
(2.2) 

for some constant c > 0. 

Consider the operator 5*^ = 5^" S^"^-^ ■ ■ ■ S^^ . Evidently, 5^ is compact, as a product of compact 
operators. To study the spectral properties of , we make use of the infinite-dimensional Perron- 
Probenius theory, as presented in [KR]. 

Lemma 1. For any f e PC[0, 1], 

(2.3) 5^/ = fl^lKf)v:^ + Af, 

where /j, > 0, tp is a strictly positive functional, S^u^ = nu^ , = 1, Au^ = 0, il>{Af) = and 

lim„_^P"||i/« <^. 

Proof. Recall that we have assumed without loss of generality that Pm = 0- Let s be the largest index 
satisfying relations Pm+i-s' = for 1 ^ s' ^ s. Define a linear operator R: PC[Q, 1] — > PC[0, 1] by 
Rf{x) = (1 — xY f{x), and let TZ be the image of PC[0, 1] under R. We consider 7?. as a Banach 
space in the induced norm, that is, ||/||7^ = maxo^x^i (f^-x)-' • 

Let K be the cone of nonnegative functions in Ti, that is, / e 7?. belongs to K \i f{x) ^0 for 
x e [0,1]. Evidently, the closure of the linear span of K coincides with TZ. We say that a linear 
operator A: TZ ^ TZ is strongly positive if it preserves K and takes any boundary point of K to an 
inner point. It is easy to see that / € is an inner point of K if liminfa,^! f{x) /(I — a;)* > 0. 

Consider the restriction of to IZ, which we denote by S^. Evidently, is compact and 
strongly positive. Hence, by Theorem 6.3 of [KR], the maximal eigenvalue fi of Sl^ is positive, 
the corresponding eigenspace is one-dimensional, and all the other eigenvalues lie strictly inside the 
circle of radius ^,. Since takes the whole PC[0, 1] to TZ, the same ^ is the maximal eigenvalue of 
S'^ ^ and all the above properties of 5^ remain valid for . It follows from the proof of Theorem 6.3 
in [KR] that (2.3) holds. □ 

By differentiating equation S^u^ = fiu^ m times we see that satisfies equation iiu'^"^^ = 
{—lyu with the boundary conditions u^'^(O) = if Pm-i = and u^'^(l) = if Pm-i = 1- Put 
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\ = II 1/™; then A is the minimal positive solution of the equation det M^(^f) = with = {—ly 
and 



^ ^m,0 ^m,m— v-i ^m,m—ir-2 

\ ^Tn,ii ^Tn,m+ii—ir-i ^Tn,m+ii—ir-2 



y^m^ir— i—ii 
^'m,ir-2—ii 

Vm,0 ) 



with ^m.j = fm.j{t)- It is easy to see that can be obtained from by the transformation 
{k, I) (r + 2 — k,r + 2 — I), and hence A is the minimum absolute value of the solutions of the 
equation detM^(t) = 0. 

Observe that the condition on the norms of ||A"|| in (2.3) means that there exists k such that 
< ij.''. Denote a = < 1 and T = (iS'^)^ 

We are approximating the operator with the help of operators defined by 



^n+l — 

It follows immediately from (2.2) that 
(2.4) 



qPm qPm-1 qPl 

'^(n+l)m'^(n+l)m-l ' ' ' '^nm+1- 



Finally, to approximate T define operators T„+i = ^5'(^+i)fc'S'(^+i)fe_i 
the above inequality that 



■S'l^ii; it follows from 



\Tn-T\\^- + o' ^ 



n 

For any function / define a sequence {/„} by /„ =Tnfn-i with /q = /. 

Lemma 2. Let {||/n||} be bounded, then the sequence fn = /n/||/n|| converges to as oo. 

Proof. By (2.3), each of /, can be uniquely represented as fi = /" + ff^, where /" is a multiple of 
and belongs to the maximal subspace invariant under A and not containing u^. Therefore, 

fn = Tnfn-1 = (T + (T„ - T))(/^i + f^l) = /^l + ^V^-l + {Tn ' T){fl_, + ftl), 

which together with = 1 gives 

ll/„"ll^ll/:-ill-||(T„-r)/„_i||, 

WfnW < + ll/n"-lll • \\{Tn " T)ti^|| + Wn " T)/^-l||. 

Recall that {||/n||} is bounded, and hence 



lib 



where /3 < 1 can be chosen arbitrary close to 1 for n big enough. Therefore, 



Wf^W < ll/n-lll fc" 



\\m\ m:-i\\ \ n 



/3| 1/^-1 1 1 



c" I a c" 
+ \^ + 



(3n 



a /3n 



5n-l + O 



Therefore, either 5n-i ^ n or (5„_i > n and 



5n < 
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for some constant a' < 1. In any case S,, ^ as n — > oo, and hence the sequence 
converges to a multiple of u^, which implies the convergence of to . □ 

For an arbitrary finite sequence a = {ai, . . . , Uk} define a piecewise constant function d G PC[0, 1] 
whose value equals ai on the interval [^^, |)forz = l,...,fc— 1 and on [^^, l] ■ Let jjfe = 
{ttfe.i, • • • , Sfe.fe}- It is easy to see that j|fe+i = kS^'"ik, which means that 



1+1 



Observe that the sequence of bmctions g„ = ('^„-)i is boimdcd. Indeed, ||jjmn 

The exponential generating function for the numbers tlmn+i,i is calculated in Corollary 1. Since the 
numerator of the corresponding expression is a polynomial in Olivier functions, which converge in 
the whole plain, the numbers ^'^^^j^ grow asymptotically as , hence for n big enough one has 

(mn)! ^ t'' ^^'^ therefore the sequence {.g„} is bounded. Moreover, (2.5) can be rewritten 

as X"^ S^_^_-^gmn+i = 5m(n+i)+i- Therefore, Lemma 2 applies, and gnm+i ^ u'^ as n ^ oo. 
Combining the above results we get 

aP \mn+l 
1- Hmn+1,1 

^™ — — = c„,iwrw, 

n^oo, ^^^i -^t (mn): 

with Cm.i = 7', and hence Theorem 2 is proved for I = 1. 

To get the proof for the other values of I one has to consider, instead of 5^, a different operator: 
gpm+i-i gPm+i-2 Its properties are identical to those of S^; to prove this one has to use 

operators R and L: f{x) 1— > x^f{x), depending on the value of Pm+i-i- D 
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